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We discuss Dyson's argument that the vacuum is unstable under a change g 2 — ► — g 2 , in the 
context of lattice gauge theory. For compact gauge groups, the partition function is well defined 
at negative g 2 , but the average plaquette P has a discontinuity when g 2 changes sign. This reflects 
a change of vacuum rather than a loss of vacuum. In addition, P has poles in the complex g 2 plane, 
located at the complex zeros of the partition function (Fisher's zeros). We discuss the relevance 
of these singularities for lattice perturbation theory. We present new methods to locate Fisher's 
zeros using numerical values for the density of state in SU (2) and U(\) pure gauge theory. We 
briefly discuss similar issues for 0(N) nonlinear sigma models where the local integrals are also 
over compact spaces. 
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1. Introduction 

Dyson instability [J], ^] - the catastrophe happening when you change the sign of e 2 in QED 
- is often invoked to limit the validity of perturbation theory and justify the factorial growth of 
the perturbative coefficients. In the functional integral formulation of scalar models, this type of 
instability is related to large field configurations 

For lattice models with compact field integration (nonlinear sigma models over compact mani- 
folds and lattice gauge theories (LGT) with compact groups), the large field problem is in principle 
absent. For g 2 < 0, the partition function is well defined and the change of sign of g 2 appears as a 
mere change in vacuum rather than a catastrophic instability. Can this explain the apparent power 
growth (rather than a factorial growth) observed in perturbative series for the average plaquette P 
in Refs. ^, [7| ||]? These series are consistent with the existence of Fisher's zeros (zeros of the 
partition function in the complex coupling plane) close to the real axis [0]. It seems clear that a 
complete knowledge of the location of the Fisher's zero would provide a compete understanding 
of the complex singularities of P. The volume dependence of these zeros also provides important 
information regarding the order of possible transitions or the absence thereof [[K], 11]. 



In these proceedings, we report recent results concerning these questions. Dyson instability is 
reviewed in Sec. ||| in the context of LGT. New "topological" methods to locate Fisher's zeros in 
SU (2) and U(l) using numerical calculations of the density of states Jl2| , |l4| ] are discussed in Sec. 
|H Similar questions for 0(N) sigma models in the complex t' Hooft coupling plane are briefly 
discussed in Sec. 0. Details can be found in a recent preprint [15]. 



2. Dyson's instability versus compact integration 

Dyson's argument goes as follows [|], §]. Suppose that a physical quantity in QED can be 
calculated as a perturbative series F(e ) = a$ + a\e 2 + . . . . If we assume that the series has a finite 
radius of convergence, then, for e 2 sufficiently small, we can interpret F{— \e 2 \) as the value of this 
quantity in a fictitious world where same charge particles attract. But in this fictitious world, every 
physical state is unstable. So, the radius of convergence is zero. Quoting the author "The argument 
[...] is lacking in mathematical rigor and in physical precision. It is intended to be suggestive, to 
serve as a basis for further discussions". 

The connection between asymptotic series and the problem of integrating large fields contri- 
butions can be understood with this very simple example 

rd<t>e-^-^ + £ t*£ / + °V-^V« . (2.1) 

The sum and the integration have been interchanged illegally. The peak of the integrand of the q-th 
order term of the r.h.s is reached when 2 = 4q. The approximation of e~ XlSf by an expansion of 
order q in A0 4 is good provided that X<p 4 «q, but at the peak of the integrand, 4 = I6q 2 and we 
need X I6q 2 «q, which fails for q large enough. On the other hand, if we introduce a field cutoff, 
as the order increases, at some order, the peak moves outside of the integration range and there 
is no factorial growth. The general expectation is that for a finite lattice, the partition function Z 
calculated with a field cutoff is convergent and ln(Z) has a finite radius of convergence controlled 
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by the zeros of the partition function. The field cutoff <j) max is an optimization parameter fixed using 
strong coupling [Q], for instance. 

A fact that is obvious but which importance regarding weak coupling expansions may have 
been overlooked is that lattice gauge theories with a compact group and nonlinear 0(N) sigma 
models have a build-in large field cutoff. In lattice gauge theory, the group elements associated 
with the links are integrated with dUi the compact Haar measure. Our notations are as follows: 
N c is the number of colors, S = Y*piaq.{^ ~ (1 /N c )ReTr(U p )) and j8 = 2N c /g ■ The number of 
plaquettes is denoted jY p = L D D(D — l)/2 . The average plaquette: P(j8) = {l/J/p) (S) will be 
our main object of study. The partition function Z(j3) is the Laplace transform of n(S), the density 
of states: 

Z(fi)= / n(S) e^ 5 , (2.2) 

Jo 

with 

n(S)=H j dUiS(S-'£(l-(l/N c )ReTr(U p ))) . (2.3) 

Assuming that ln(n(5)) is extensive we can write 

n{S)=t^ s '^ . (2.4) 

It is important to notice that at finite volume, S max is finite. For instance, S max = 2jV p for SU (2N) 
and \ jV p for SU (3). In the strong coupling expansion, we expand in power of /3: Z = L~ = oZ«j3" 
with \z n \ < S'^ax/nl, so at finite volume, Z is an analytical function, not only on the negative real 
axis, but over the entire j8 plane. 



On the other hand, it is possible to show that for SU {2N) on even lattices Jig ] 

Z(- j 8)=e 2 ' 3 ^Z( j 8) . (2.5) 

Consequently, 

n(2 r yy p - S) = n(S) and P(j8) +P(-j8) = 2 (2.6) 

Since limp^ +00 P(j3) = 0, P has a discontinuity at g 2 = and a regular series for P about g 2 = is 
not possible. However, it does not necessarily mean that the series has factorial growth. 



It is useful to consider first the case of a single SU(2) plaquette [17]. In that case, n(S) = 
^ \/S(2 — S) (invariant under S ^2 — S). The large order of the weak coupling expansion j3 — > +°° 
is determined by the behavior of n(S) near S = 2, itself probed when j3 —> — °° in agreement with 
the common wisdom that the large order behavior of weak coupling series can be understood in 
terms of the behavior at small negative coupling, y/2 — S is then expended about S = (radius of 
convergence = 2). This yields the convergent expansion 



/to /!(l/2-/)7o 

As expected this is a not a regular series in the sense that the "coefficients" of /3~' depend on j3, but 
in a way that is invisible in perturbation theory. The crucial step is to get j3 -independent coefficients 
by neglecting the missing tails of integration. 

p2B poo 

/ dtz-'t'+'l 2 - / Ae-y+ 1 /2 + <9( e - 2 ' 3 ) (2.8) 
Jo Jo 
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which in turn creates a factorial growth of the coefficients. The peak of the integrand crosses the 
boundary near order 2/3. Dropping higher order terms (than order ~ 2/3) agrees with the rule of 
thumb (minimizing the first contribution dropped). The non-perturbative part can be fully recon- 



structed (higher orders + tails) [18]. For L 4 lattices, the crossing should be near order 2fijVp. 
Non-perturbative effects should be explainable by the contributions near S max . We plan to study 
this question on small lattices. 

3. Fisher's zeros from the density of states 

The poles of P are located at the Fisher's zeros. At finite volume, we expect these zeros to be 
isolated in the /3 plane. It seems plausible that the zeros will accumulate along lines going through 
in the 1 / jS plane as they do for Bessel functions. It is possible to use n(S) to calculate Z at complex 



/3. The calculation of n(S) for SU (2) is discussed in Ref. [12]. Additional checks were made by 
calculating the first three moments of n(S). For U(l) lattice gauge theory, multicanonical methods 
relying on the Biased Metropolis-Heatbath Algorithm Jl3| ] were used [[l4|]. Using this U(l) density 
of states, we have calculated the plaquette distribution calculated at fixed /3 and checked that there 
is an approximately symmetric double peak near B = 0.979 for a 4 4 lattice. 

For both U(l) and SU (2) on a 4 4 lattice, the numerical calculation of Z(/3) with Imfi ~ 0.2 is 
difficult because /3 is multiplied by ,jV p and the integrand oscillate rapidly. A preliminary idea of 
the distribution of zero can be obtained using semi-classical methods. Using the "color entropy" 



f(s) defined in Eq. (2.4), the saddle point of the integral is at sq given by solving f'(so) = j3. 



Z becomes a Gaussian integral with correction of order yT/ ' jV p as long as Ref"(so) < 0. As a 



Gaussian density of states has no complex zeros [19], it seems clear that zeros should appear in 
regions of the /3 plane corresponding to regions of the s plane such that Ref'(so) > 0. Using 
Chebyshev approximations of f(s), we have constructed the boundary (Ref"(s) = 0). The results 
are shown in Fig. |[ The boundary form narrow tongues ending at a complex zeros of /". These 
complex zeros are then mapped in the /3 plane using /'. Their number depends on the degree of the 
polynomial approximation, but the general shape is robust under changes in the degree. It appears 
that in the case of SU (2) the images in the /3 plane are never on the real axis in contrast to the case 
of 1/(1). 

New methods have been developed to locate the Fisher zeros [15, ^C| ]. Given the fact that Z is 
an entire function in the /3 plane, and that P = — (dZ/dp)/Z, the worse thing that can happen to P 
is that Z has a zero of order k, say at /3o- Then (dZ/dfi)/Z ~ k/(fi — /3o) for /3 ~ /3o- If we now 
integrate over a closed contour C, 

{iln)- 1 I dp{dZ/dp)/Z = Y n k(C) , (3.1) 
Jc T 

where n^C) is the number of zeros of order k inside C . This allows us to monitor the accuracy of 
the calculation. We need to check that in good approximation, the real part is an integer and the 
imaginary part is zero. This is illustrated in Fig. |2] for a rectangular contour of variable height in 
the p plane. Despite these encouraging results, there remain dependence on the interpolation or fit 
used to evaluate f(s) numerically. Resolving this issue should allow us to find finite size scaling 



for the zeros as discussed in Refs. [10, 11]. 
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zeros of f" (SU(2)) zeros of f" (U(l)) 
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f ' (s) for zeros (SU(2)) f'(s) for zeros (U(l)) 
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Figure 1: Top: complex zeros and zeros of the real part of f"{s) in the complex s plane with 40 Chebyshev 
polynomials on 4 4 for SU (2) (left) and U(l) (right). Bottom: f'(s) evaluated at the complex zeros of f"(s) 
shown on the previous figure for SU(2) (left) and t/(l) (right). 
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Figure 2: Re (left) and Im (right) part of Y,k n k defined in Eq. (3.1) for a rectangular contour with 2.1 < 
Refi < 2.3 and < /mj8 < y with a variable y, for SU(2) on a 4 4 lattice. Two independent numerical values 
of n(S) were used. 
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Figure 3: Left: images in the A' plane of lines of constant imaginary part 2.25, 1.75, 1.25,0.75,0.25,-0.25, 
....,-2.25 in the complex mass gap plane and of the singular points (red dots) for a 8x8 lattice. Right: Fisher 
zeros for N — 2 (blue) and images of singular points (red). 

4. 2 — D 0(N) nonlinear sigma models 

The nonlinear 0(N) sigma models on even cubic lattices have similar properties under the 
exchange of the sign of the coupling, namely Z{— g 2 ] = e 4DL °/ g2 Z[g 2 ]. The complex singularities 
of the average energy in 2 dimensions, for complex 't Hooft coupling X' = g^N have been studied 
in the large-A^ limit . Details can be found in a recent publication []TJ]. A striking difference with 
the linear model is the absence of cut along the negative real axis. It was argued that the Fisher's 
zeros can only be inside a clover shaped region of the complex V plane or equivalently outside of a 
region delimited by 4 approximate hyperbolas with asymptotes on the boundary of a cross of width 
0.5 centered at the origin in the 1 /X* plane. The argument holds for large N and large volume. This 
limit is being studied using exact results at finite ,/V and finite volume. The graphs of Fig. || made 
with N = 2 on a 8 2 lattice can be compared with the corresponding ones in Ref. [|l5|]. 

5. Conclusions 

For pure gauge models with compact groups, there is no loss of vacuum when g^ — ► — gg, but 
only a change of vacuum. The discontinuity of the plaquette forbids the existence of a converging 
perturbative series but does not dictate the large order behavior. Reliable methods to locate Fisher's 
zeros are in progress. Non-perturbative effects should be accountable by modified expansions. New 
data for perturbative coefficients should help in this task . 
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